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Fry (1996) showed that galaxy bias has the tendency to evolve towards unity, i.e. in the long 
run, the galaxy distribution tends to trace that of matter. Generalizing slightly Fry's reasoning, 
we show that his conclusion remains valid in theories of modified gravity (or equivalently, complex 
clustered dark energy). This is not surprising: as long as both galaxies and matter are subject to 
the same force, dynamics would drive them towards tracing each other. This holds, for instance, in 
theories where both galaxies and matter move on geodesies. This relaxation of bias towards unity is 
tempered by cosmic acceleration, however: the bias tends towards unity but does not quite make it, 
unless the formation bias were close to unity. Our argument is extended in a straightforward manner 
to the case of a stochastic or nonlinear bias. An important corollary is that dynamical evolution 
could imprint a scale dependence on the large scale galaxy bias. This is especially pronounced 
if non-standard gravity introduces new scales to the problem: the bias at different scales relaxes 
at different rates, the larger scales generally more slowly and retaining a longer memory of the 
initial bias. A consistency test of the current (general relativity + uniform dark energy) paradigm 
is therefore to look for departure from a scale independent bias on large scales. A simple way is 
to measure the relative bias of different populations of galaxies which are at different stages of 
bias relaxation. Lastly, we comment on the possibility of directly testing the Poisson equation on 
cosmological scales, as opposed to indirectly through the growth factor. 

PACS numbers: 98.80.-k; 98.80.Bp; 98.80.Es; 98.65.Dx; 95.35.+d; 95.36.+x; 95.30.Sf 



I. INTRODUCTION 

In a classic paper, Fry [l[ showed that galaxies, once 
formed, would evolve to eventually trace mass. This re- 
sult was subsequently extended to the case of a stochas- 
tic bias by Tegmark and Peebles [3. Similar results 
were found in numerical simulations The linear bias 
paradigm has been applied to a variety of cosmological 
models, including ACDM and quintessence Both [l| 
and 0] assumed Newtonian/Einsteinian gravity from the 
start. One of our principal aims is to show that relaxing 
this assumption does not alter the basic conclusions. In 
the process, we will uncover an interesting, unanticipated 
corollary. 

The basic reasoning is simple: both galaxies and (dark) 
matter can be regarded as test particles in some gravita- 
tional field; since their dynamics are identical, they are 
expected to trace each other given sufficient time. This 
holds regardless of the origin of the gravitational field. 
Departures from Newtonian/Einsteinian gravity on large 
scales have been suggested as a possible explanation for 
the observed cosmic acceleration. Most of these theories 
of modified gravity, for instance DGP [5| or scalar-tensor 
theories @, are metric theories, and matter (or galaxy) 
particles move on geodesies defined by the metric. This 
is true as long as matter and galaxies are minimally cou- 
pled to gravity in the relevant frame 7|. The statement 
that matter and galaxies should evolve to trace each other 
does not rely on how the metric comes about. All that 
matters is that they have the same dynamics under the 
metric. More generally, as we will see, our conclusions on 



bias evolution remain valid as long as both galaxies and 
matter are subject to the same force. 

Note that the same conclusion holds if, instead of mod- 
ifying gravity, one introduces additional sources of fluctu- 
ations, for instance clustered dark energy. The presence 
of dark energy perturbations modifies the relation be- 
tween matter perturbations and the metric, or the grav- 
itational potential (just as modifying gravity does). But 
once this gravitational potential is given, both galaxies 
and matter particles respond to it in the same way, there- 
fore guaranteeing eventual evolution towards a trivial bi- 
asing relation between them. This is consistent with the 
fact that modifying gravity can be viewed as some com- 
plicated form of dark energy - no surprise given that one 
can always shift terms from the left hand side of the Ein- 
stein equation to the right [8|. Modifications of gravity 
are indistinguishable from complex dark energy, at least 
at the linear level, unless constraints are put on the na- 
ture of dark energy, such as spatial homogeneity or zero 
coupling to matter (e.g. (9l. [lfj] and references therein). 

Why are we interested in proving these general but in 
a sense rather obvious statements about bias evolution? 
First, as we will see, there is an interesting and proba- 
bly not so obvious corollary: that bias evolution can im- 
print a scale dependence even if the initial bias were scale 
independent. The effect is especially large if one's non- 
standard theory of gravity/dark energy introduces new 
scales to the problem, making the growth factor scale 
dependent. A promising way to test the current (general 
relativity + uniform dark energy) paradigm is therefore 
to look for a scale dependent bias. Our second (and orig- 
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inal) motivation is to address the issue of directly testing 
the Poisson equation. Most observations, such as grav- 
itational lensing and the microwave background, tell us 
the metric fluctuations or gravitational potential, the left 
hand side of the Poisson equation. There are very few ob- 
servations that tell us directly what sources the gravita- 
tional potential, the right hand side of the Poisson equa- 
tion. The galaxy distribution provides one such handle, 
but it is complicated by the uncertain galaxy bias. One 
might hope that the tendency of the galaxy bias to relax 
towards unity could help bridge the gap between galax- 
ies and mass. However, as we will see, in an accelerating 
universe, the bias relaxation process is quite slow. Modi- 
fications of gravity that predict cosmic acceleration tend 
to slow it further. The slow relaxation process makes 
it difficult to argue that galaxies that formed at high 
redshifts tend to have a trivial bias today, unless they 
formed with a bias close to unity. On the other hand, as 
we will see, the slow bias relaxation is a boon to the scale 
dependence test (motivation number one). 

Two clarifying remarks are in order. First, throughout 
this paper, the galaxy bias is defined with respect to the 
matter fluctuations. It is important to bear this in mind 
since, in models that allow clustered dark energy, there 
are different kinds of galaxy bias: with respect to mat- 
ter, dark energy or the sum. From the point of view of 
dynamical evolution, the bias with respect to matter is 
the natural one to consider. Secondly, an important as- 
sumption in [l[ as well as here is that the galaxy number 
is conserved, just as mass is. In practice, galaxies are of 
course created (they form) and destroyed (they merge). 
A complete theory of galaxy bias needs to fold in these 
processes. We will have more to say about this in fJVT] 
Dynamical evolution is in any case an important part of 
the story. 

This paper is organized as follows. The basic equations 
are reviewed in i flll The implications for the linear galaxy 
bias are worked out in ^IIII The extension to a (linear) 
stochastic bias is presented in ^IVl and the extension to 
a nonlinear bias is given in ifVl The interesting corollary 
of a possible scale dependent imprint on the galaxy bias 
is derived and discussed at the end of both ^IIII and IIVI 
In £|VI[ we conclude with some remarks about testing 
theories of modified gravity or complex dark energy. 

II. BASICS 

According to the current structure formation 
paradigm, the fundamental equations governing the 
dynamics of sub-Hubble matter fluctuations are [l2j: 

S' m = -V-[(1 + S m )v m ] (1) 

Vm + (Vm ' V) V m + -V m = -V(j> (2) 

a 

V 2 = Ai:Gp m a 2 5 m (3) 

where 5 m = (p m — p m )/Pm is the matter overdensity, with 
p m being the matter mass density and p m its spatial av- 



erage, v m is the matter peculiar velocity, a is the scale 
factor, <f> is the gravitational potential or metric perturba- 
tion, G is Newton's constant, ' denotes a derivative with 
respect to conformal time 77, and V denotes a derivative 
with respect to comoving position x. 

Equations (fTJ) and ^ express respectively mass and 
momentum conservation. Equation ([3]) is the Pois- 
son equation, which tells us how matter fluctuations 
source the gravitational potential according to Newto- 
nian/Einstcinian gravity. 

Galaxies (or more precisely, some population thereof), 
once formed, obey the following equations on sub-Hubble 
scales: 

6' g = -V-[(l + 5 g )v s ] (4) 
a' 

v g + (v g • V) v g + -v g = -V0 (5) 

where S g = (n g — n g )/n g is the galaxy overdensity, with 
n g being the galaxy number density and n g its spatial 
average, and v g is the galaxy peculiar velocity. Equation 
(j4]) crucially assumes number conservation, a subject to 
which we will return in WII Equation (J5]) is momentum 
conservation, or the equation of motion, for the galaxies. 
It is worth noting that we do not assume v m = v g , as 
is often done. This is in part because we want to make 
explicit the importance of subjecting the galaxies and 
matter to the same force. If they are not subject to the 
same force, v m and v g would diverge from each other; if 
they are, the two velocities tend to converge, as we will 
see. 

Equations |T]) to ([5]) constitute the starting point for 
[HQ on the evolution of bias. Our key observation is that 
Eq. §5§ is unnecessary, at least on a qualitative level. 

This is interesting in light of the recent interest in 
theories of modified gravity. By and large, these theo- 
ries respect mass and momentum conservation; particles 
(matter or galaxy) move on geodesies. The only equation 
that is modified in these theories is therefore the Poisson 
equation ([3]). Similarly, in theories where dark energy 
has non-trivial clustering behavior, dark energy can act 
as an additional source for gravitational perturbations, 
and therefore the simple form of the Poisson equation 
given in ([3J is again modified. 

Henceforth, our derivation relies on Eqs. (H}, Q 
and ([5]), but not on (J3J). In other words, our only assump- 
tions are: that mass and number are conserved, and that 
matter and galaxies are subject to the same force. 

III. LINEAR BIAS 

Keeping only terms linear in fluctuations, one can see 
from Eqs. © and © that 

(v g - v m )' + — (v g - v m ) = (6) 
which implies 

(v g - v m ) = u /a (7) 
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where Uo is independent of time but dependent on posi- 
tion. Note that in deriving this result, we do not need the 
right hand sides of Eqs. ((2]) and (0 to be the gradient of 
some potential; all we need is that they are identical i.e. 
matter and galaxies are subject to the same force. Equa- 
tion ([7]) tells us that if v g = v m initially, it will remain 
so; any initial velocity bias, should it exist, decays away. 
This holds in the context of linear perturbation theory 
(i.e. large scales) [Tlj |. In fact, [l|, I2J set the galaxy and 
matter velocities to be equal from the outset. 

Similarly, linearizing Eqs. ([TJ) and ^ and taking the 
difference, we obtain 

(S g - 5 m y = -V • (v g - v m ) = -V • u /a . (8) 

This can be integrated to give 

f da 

5 g - S m = -V • u / + A (9) 

where Ao is independent of time but spatially dependent, 
and H is the Hubble parameter. The integral over a gen- 
erally gives a decaying term in contrast with the constant 
A term, and so to leading order, the galaxy bias evolves 
as: 

o m (a) o m (a) 

where we have inserted the argument a in each quantity 
that is time dependent (while space dependence is kept 
implicit). The quantity Ao is independent of time, and 
can be alternatively expressed as 

A = {b - l)S m (a ) (11) 

where bo is the initial bias at some early time correspond- 
ing to scale factor ao- From Eq. (fTT))) . we can see that as 
long as the matter fluctuation S m grows, the galaxy bias 
b evolves towards unity. Exactly how rapidly it evolves 
depends on the growth rate of S m , which cannot be deter- 
mined unless the modification to Eq. ([3]), or lack thereof, 
is specified. For small fluctuations, it is not unreason- 
able to expect that S m grows roughly on the Hubble time 
scale. 

There is an interesting corollary that follows from Eq. 
(jlOp . Suppose the galaxy bias is initially scale indepen- 
dent. This means <5 s (ao,x) = 6o<5 m (ao,x), where bo is 
the initial bias which is independent of position, and do 
denotes the scale factor at some early time. Equations 
(fTT))) and Hi]) then tell us for a > a : 

hi \ -i , / 7 -,^m(ao,x) 

&(a,x) = l + (&o-lK— ; r 12 

d m (a,x) 

which means b is independent of position, or is scale in- 
dependent, only if (5 m (a,x) preserves exactly the initial 
spatial dependence encoded in <5 m (ao,x). We know this 
holds in general relativity (plus uniform dark energy) at 
the level of linear theory in the sub-Hubble limit. One 



can define the bias b in Fourier space as well, in which 
case Eq. (fT2"j) takes the form: 

6(a,k) = l + (fe -l)§4 £ #- ( 13 ) 
d m (a,k) 

In the standard paradigm of general relativity plus uni- 
form dark energy, different sub-Hubble Fourier modes of 
8 m grow with the same (k independent) growth factor, 
hence an initially scale independent bias (bo) will remain 
scale independent (b). Note that according to the stan- 
dard paradigm, the initial linear bias is expected to be 
scale independent [H, 0, [H| . 

However, it is entirely possible that modifications to 
the Poisson equation ([3]) take a scale dependent form, 
which would make the perturbation growth scale depen- 
dent. General relativity itself introduces such correc- 
tions, but they are generally suppressed by (aH/k) 2 and 
are therefore very small except close to the Hubble scale 
where cosmic variance is large [ill ]. A less trivial exam- 
ple is a Yukawa-like modification of the Poisson equa- 
tion: V 2 — > V 2 — /i 2 , where 1/fi defines some new scale 
[l7l EH- Scalar-tensor theories introduce modifications of 
this kind where fj, 2 is related to the second derivative of 
the scalar potential [l9j , in addition to the better known 
effect of giving G a time dependence [2(| H2, H3|- ^ n 
these cases, the linear growth factor S m (a, k)/5 m (ao, k) 
will be k dependent. As a concrete example, we show 
in the bottom two panels of Fig. [JJ the linear bias as a 
function of redshift and wavenumber k for two different 
populations of galaxies, with the following modification 
to the Poisson equation: 

(V 2 - a 2 m 2 )<t) = 4irGpa 2 S m (14) 

where to is a constant and a 2 m 2 plays the role of /i 2 . 
We use fj, oc a such that the Yukawa correction becomes 
small in the early universe [l8| , although some other time 
dependence is possible. There is also the question of 
how such a scalar-tensor theory can be made consistent 
with solar system constraints; this depends on the full 
Lagrangian of the theory - the kinetic term, the poten- 
tial and the coupling to curvature can always be cho- 
sen to satisfy such constraints [23|. In Fig. [TJ we adopt 
to = 0.05 h/Mpc which appears to be consistent with cur- 
rent observations A flat cosmology of Q m = 0.27 
and £1a = 0.73 is assumed in computing the background 
expansion rate. 

From Eq. (fT4|) . one can see that the high k modes 
(k/a 3> to) are insensitive to the Yukawa correction. 
Hence, the high k end of Fig. [JJ shows a bias evolution 
that is consistent with Einsteinian/Newtonian gravity. It 
is worth stressing that while the bias wants to evolve to- 
wards unity, it does not quite make it to unity. This is 
the result of cosmic acceleration: perturbation growth 
is much suppressed once the universe starts to acceler- 
ate, and the bias relaxation process is stalled. In fact, 
even if we continue the evolution into the future, the 
bias will not become very close to unity, unless its ini- 
tial value were close to unity to begin with. According 
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FIG. 1: The evolution of the linear bias and relative bias 
as a function of wavenumber k. In the bottom panel is b A , 
the linear bias of galaxies that formed at redshift 5 with an 
initial scale independent bias of bo ~ 2. In the middle panel is 
b B , the linear bias of galaxies that formed at redshift 2 with 
an initial scale independent bias of bo = 0.7. The top panel 
shows the ratio of the two. The number labeling each curve 
is the corresponding redshift. This makes use of the Yukawa 
modification to the Poisson equation: Eq. (|14|l . 



to general relativity, S m (a = oo)/5 m (a = 1) ~ 1.36 for 
our adopted parameters, implying bias relaxation via Eq. 
(I13p is rather limited in the future. 

If general relativity were the whole story, the curves 
in Fig. Q] would have been flat (aside from very small 
(aH/k) 2 corrections mentioned earlier). Instead, we see 
that the Yukawa correction further suppresses the growth 
of perturbations at low /c's (compared to high fc's), which 
significantly slows their evolution towards a bias of unity. 
This is qualitatively what one expects from theories of 
modified gravity which are constructed to yield cosmic 
acceleration: such theories tend to make gravity weaker 
on large scales. In other words, the lower-A: modes retain 
a longer memory of the initial bias, while the higher- 
k modes catch up more efficiently (but not completely) 
with the matter perturbations. The result is a scale de- 
pendent bias, as shown in Fig. [T] Note that, in prin- 
ciple, this scale dependence is temporary - if the bias 
does eventually reach unity in the far future, it becomes 
scale independent once again. In practice, in an acceler- 
ating universe, the bias relaxation process is so slow that 
this does not actually occur, unless the initial bias were 
quite close to unity. This is a useful feature from the 
perspective of preserving, and potentially observing, the 
dynamically generated scale dependent bias. 



FIG. 2: Analog of Fig. [T]except that a DGP motivated mod- 
ification of growth rate is used: Eq. (|15[) . Note that this does 
not account for the scale dependence arising from the tran- 
sition from the scalar-tensor regime to the general relativity 
regime, the so called r, effect. 



We note in passing that the rather popular f(R) the- 
ory is a scalar-tensor theory, and is therefore expected 
to exhibit a scale dependent bias qualitatively similar to 
that shown in Fig. [1] It appears, however, that the par- 
ticular f(R) theories usually considered in the literature 
have an effective mass m (or /i) that is rather small [||, 
which makes for a weak scale dependence. There is no 
fundamental reason, though, to focus only on these par- 
ticular realizations of a scalar-tensor theory. 

Another example of a modified gravity theory is DGP 
H, which is expected to introduce scale dependence to 
the growth factor as well. There are two sources of scale 
dependence. One is that DGP is expected to give correc- 
tions to the growth factor of the order of aHo/k, arising 
from its peculiar nonlocal character from the 4D perspec- 
tive (Scoccimarro, private communication). The other is 
that DGP introduces a special scale, often referred to 
as r*, where gravity makes a transition from the gen- 
eral relativity regime to a scalar-tensor regime 24] . One 
expects the structure growth to have a feature around 
r*. Its precise value is however uncertain, because it de- 
pends on the exact nature and characteristics of the fluc- 
tuations. None of the existing calculations of structure 
growth in DGP gravity address growth in this transition 
regime |25j . The general expectation is that is close 
to the nonlinear scale, but its actual value could well be 
larger. As an illustration, in Fig. [2J we sidestep the r* is- 
sue, and simply modulate the DGP growth by an aHg/k 
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correction (see also jig]): 



IV. (LINEAR) STOCHASTIC BIAS 



5%*- (a, k) = C sson («, k)(l - aHo/k) (15) 

where S^ lsson denotes the matter fluctuation according 
to the Poisson equation, with a background expansion 
chosen to match the expected DGP growth on small 
scales (but larger than r*) [2(|. Here, the scale depen- 
dence of the bias is much weaker than in the Yukawa 
example, though one must remember that this ignores 
the r* effect. As mentioned before, general relativity it- 
self introduces corrections of this sort, but they come in 
at quadratic order (aH/k) 2 and are even weaker. 

Lastly, instead of modifying the Poisson equation ([3]) 
on the left hand side, one could also modify it on the 
right hand side by having dark energy that clusters. This 
would introduce new scales to the problem, such as the 
Jeans scale, and one expects the growth of perturbations, 
and therefore the bias relaxation, to become scale depen- 
dent. 

We should add that in theories where gravity is mod- 
ified, or clustered dark energy is introduced, it is quite 
possible that the galaxies are born with a scale dependent 
bias on large scales, contrary to expectations according 
to the standard excursion set theory of halos [H, 03] . It 
would be very surprising if any such initial scale depen- 
dence conspires to cancel the scale dependence that de- 
velops dynamically. Checking for departure from a scale 
independent bias on large scales therefore constitutes an 
interesting test of (scale dependent) modifications to the 
Poisson equation (J3]) . One way to perform this test is to 
compare the galaxy power spectrum and the mass power 
spectrum from lensing (the latter requires the Poisson 
equation to convert from metric to mass fluctuations) . A 
simpler test is to compare the large scale bias of different 
populations of galaxies, which are at different stages of 
dynamical bias relaxation. For instance, labeling the two 
populations of galaxies by A and B, their relative bias is 



b A (a,k) _ l + (6^-l)5 m (a^k)/5 m (a,k) 



b B (a, k) 



(&£-l)Ma£,k)/<Wa,k) 



(16) 



for a > ciq. A ratio of this kind can be obtained from the 
ratio of the observed power spectra of A and B. The top 
panels of Fig. [1] and [2] illustrate this ratio for two par- 
ticular populations of galaxies. It is worth noting that 
the best chance for observing a scale dependent relative 
bias is to compare overbias (> 1) galaxies with underbias 
(< 1) galaxies. This is because modified gravity theories 
tend to predict a higher bias for the overbiased galaxies, 
and a lower bias for the underbiased galaxies. This effect 
gets stronger at lower fc's; the scale dependence is most 
visible in their relative bias. So far, there is no observa- 
tional evidence for departure from scale independence on 
large (i.e. linear) scales [271 ] . 



Having a stochastic bias means one introduces an ad- 
ditional parameter into the problem, following Dekel & 
Lahav [28|]. In addition to b, defined by, 



b 2 = {S 2 g )/{S 2 J 
one introduces r, defined by 

br ee (S g 5 m )/(5 2 m ) 



(17) 



(18) 



Here, one can think of the ensemble averages in several 
different ways. For instance, (S g S m ) can be thought of 
as the ensemble average of the product of 5 g and 5 m , 
each smoothed on some given scale R. It can also be 
thought of as the two point correlation function of 5 g 
and S m at some given comoving separation Ax. Another 
alternative is to think of it as the power spectrum at a 
wavenumber k: oc [(5 g (k)S^(k)} + (<5*(k)6 m (k)}]/2. De- 
pending therefore on the definition, b and r are in general 
functions of R, Ax or k respectively. It is worth noting 
that cross-correlation measurements of galaxies suggest 
r is consistent with unity on large scales [27j ■ 

Squaring Eq. ^ , and keeping only the leading contri- 
bution to the ensemble average, we have 



b 2 -2br + l = (A 2 )/(6 2 m ) 



(19) 



where (Aq) is time independent, but (S 2 ^) is not. Alter- 
natively, multiplying Eq. ^ by 8 m , and keeping only 
the leading contribution: 



&r - 1 = (A <W/<<O 



(20) 



The evolution of the linear galaxy bias b and the cross- 
correlation coefficient r is therefore given by: 



, (Ao<U , (Aj) 
(A <5„ 



1/2 



(SI) 



(21) 
(22) 



Recall that A is independent of time (Eq. [5]). As- 
suming the matter fluctuation S m grows, we can see that 
at late times: 



A2) (A 6 m ) 



(S 2 m ) (S? n ) 

and therefore at late times: 

(A <5 m ) 

(SI) 
MA§> , 



1 



1 



(A S m ) 
(SI) 



(23) 

(24) 
(25) 



and so, b — > 1 and r — * 1, as expected, with r tending 
to unity at a faster rate than b, although, as emphasized 



6 



before, cosmic acceleration slows the approach to unity 
considerably. 

Our remarks in ifllll concerning the possible develop- 
ment of a scale dependence to the large scale galaxy bias 
apply to b and r here as well. In other words, the quan- 
tities that control the evolution of b and r, (Ao8 m ) / (8^) 
and (Aq)/((5, 2 ti ), can develop a scale dependence (on R, 
Ax or k), due to possible scale dependent modifications 
of the Poisson equation. 



V. NONLINEAR BIAS 

Fry & Gaztahaga (29l | introduced a systematic expan- 
sion that relates the galaxy and matter fluctuations: 

5 g = bS m + jC + - (26) 

To find the evolution of b 2 , we need to keep terms up to 
second order in fluctuations. Assuming gradient flow i.e. 
v m = V(j)" m and v g = V<#J, Eqs. © and © tell us that 

(^-C)'+-(^-C) (27) 

+iv(^-C)-v(^ + C) = o- 

To linear order, this equation implies what we already 
know from Eq. (J7), i.e. (4> v g — 4> v m )^ decays like 1/a, 
where the superscript W denotes the order we are con- 
sidering. Whether (cj> v — decays depends on the 

behavior of {4>Z + <t>m) which depends on how the Pois- 
son equation ([3]) might be modified. It can be shown that 
(< #{ — <t >™)^ decays away if (<?^+<^m) grows slower than 
a' |3C| . Henceforth, we will assume so, following [l| who 
essentially assumed v g = v m from the beginning. 

Equation (fTJ) can be rewritten, up to second order, as 

-V-v m (x) = <4(x) (28) 

~ / (2^3 (2^ elk X "fc^ <5 " l(k),5m(k ~ k) 

where all time dependence is left implicit, and k and k 
denote wavevectors. We have abused the notation a bit 
by using 8 m to stand for both the real space and Fourier 
space counterparts, distinguishing them purely by their 
arguments. A similar equation holds for the galaxies as 
well. Assuming the difference v g — v m decays away to 
second order (or is zero to begin with), we have therefore 

few -«*>r = /J||^'^ <»> 

KWC-'l-WWi-l] 

which holds to second order. Expanding 8 g as in Eq. 
(|26[) . the first order terms in the above expression tell us 



that (b— 1) decays like 1/Sm , consistent with the finding 
in mill The second order terms go like 

l(b-l)6$ + b -l(5£>) 2 }'= (30) 

r d 3 ~k d 3 k fak-k (1) - 

J (2tt)3 (2tt)3 6 fc 2 6m W ( ' 

where we have suppressed the x dependence on the left 
hand side, and the time dependence on both sides. Here, 
Aq 1 ' is the Fourier counterpart of A in Eq. (fTO)) , and 
is therefore a first order quantity, and denoted as such. 
Note also that it is time independent. 

It is straightforward to see that, to leading order: 

b2 = 0[8<n/8$ i ] + 0[l/5g>] (31) 

One generally expects that 8 m ' ~ 8 m ' , and therefore 62 
tends towards zero at late times as long as 8^ grows. 
The same caveat applies as before: cosmic acceleration 
significantly slows the relaxation of 62- 

It is noteworthy that 62 develops a scale dependence 
even if it started off scale independent. This is true even 
if the Poisson equation © holds fj. The same is not 
true for b: the Poisson equation ([3]) guarantees b stays 
scale independent if it started off so, as noted in mill 

VI. DISCUSSION 

Our results on the evolution of bias are summarized by 
Eq. CEO]), P or for the linear bias, by Eqs. |jSI) 
and (12"2")) for the (linear) stochastic bias, and by Eq. (f3"Tj) 
for the nonlinear bias. The Poisson equation ([3]) is not 
assumed in deriving these results, therefore leaving open 
the possibility that gravity is modified or new sources of 
fluctuations (such as clustered dark energy) exist. The 
only assumptions made are mass and number conserva- 
tion, and that matter and galaxies are subject to the 
same force. For the linear bias and the stochastic bias, 
it is not even necessary that the force be the gradient of 
some potential. 

Following these assumptions, our main conclusions are: 
(1) the bias tends to relax towards unity, but this relax- 
ation is significantly slowed by the diminishing structure 
growth once cosmic acceleration commences; (2) the cor- 
relation coefficient r tends towards unity faster than the 
bias b does; (3) the large scale (linear) bias dynamically 
develops a scale dependence if modifications to the Pois- 
son equation are scale dependent. 

Conclusion (1) implies, everything being equal, galax- 
ies that formed at higher redshifts tend to have a more 
relaxed bias today. The stalling of bias relaxation by 
cosmic acceleration, however, implies we cannot auto- 
matically conclude that such galaxies have close to unity 
bias today, unless they had an initial bias that was fairly 
close to unity. In other words, one cannot freely approx- 
imate S m ~ 8 g by appealing to bias relaxation alone. 
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This is unfortunate, because S g is about the only ob- 
servational handle we have on S m without assuming the 
Poisson equation. Other observations, such as the mi- 
crowave background and gravitational lensing, generally 
probe the metric perturbations, i.e. the left hand side of 
the Poisson equation [3l|. A direct test of the Poisson 
equation requires independent knowledge of 8 m , i.e. the 
right hand side. It appears one has to be content with 
indirect tests via the linear growth factor, such as those 
widely discussed in the literature [32] . They test for how 
<f> changes with time, which the system of equations |1}, 
@ and predicts. The relation between <\> and S m is 
tested only in so far as the prediction for linear growth 
from this whole system of equations is tested. To go be- 
yond this, one would need some independent method for 
constraining the galaxy bias. Standard methods to do 
so carry additional baggage. Methods using higher mo- 
ments, for instance, assume Gaussian initial conditions 
as well as the Poisson equation itself. 

Conclusion (3) is, in our view, the most interesting one. 
It suggests a simple way to test for scale dependent mod- 
ifications to the Poisson equation, which are expected in 
many existing theories of modified gravity or clustered 
dark energy (see £ )IIII) . One can measure the relative 
large scale (linear) bias of two different populations of 
galaxies which are at different stages of bias relaxation 
(Eq. [TB] ) • The standard paradigm of general relativity 
+ uniform dark energy predicts a scale independent lin- 
ear bias (both at birth and in its subsequent evolution). 
Checking for departures from scale independence there- 
fore constitutes a test of the paradigm. As illustrated in 
the top panels of Fig. [T]and[2l this test is best carried out 
by comparing overbiased and underbiased galaxies. It is 
also important to focus on large scales, since the bias is 
expected to develop scale dependence on nonlinear scales. 

It was recently noted by [33] that primordial non- 
Gaussianity could give rise to a scale dependent large 



scale bias as well. Distinguishing this possibility from 
modifications to the Poisson equation should be feasi- 
ble since there are other more direct ways to test for 
primordial non-Gaussianity, such as from the microwave 
background anisotropies [341 ]. 

As emphasized in our calculation of the bias evo- 
lution assumes galaxy number conservation, just as in 
Ultimately, one would like to go beyond the treat- 
ment in this paper, and formulate a complete theory of 
halo bias accounting for formation, mergers and dynam- 
ics, extending the excursion set theory [l3|, [l4| beyond 
the standard (general relativity + uniform dark energy) 
paradigm. Such a theory would depend on the details of 
how the Poisson equation is modified. It is interesting 
to note that the standard excursion set theory predicts 
a bias evolution that is consistent with the one derived 
here: b(a) — 1 oc 1/ D(a) where D is the linear growth fac- 
tor (the halos are identified at some scale factor eto < a). 
This suggests our simple dynamical calculation already 
captures much of the relevant physics. An additional 
bonus of extending the excursion set theory beyond the 
standard paradigm is that one could address the issue of 
whether galaxies/halos could be born with a scale depen- 
dent bias in non-standard theories. We hope to pursue 
these issues in the future. 
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